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Future quantum technologies rely heavily on good protection of quantum entanglement against 
environment-induced decoherence. A recent study showed that an extension of Uhrig's dynamical 
decoupling (UDD) sequence can (in theory) lock an arbitrary but known two-qubit entangled state to 
the iVth order using a sequence of N control pulses [Mukhtar et al, Phys. Rev. A 81, 012331 (2010)]. 
By nesting three layers of explicitly constructed UDD sequences, here we first consider the protection 
of unknown two-qubit states as superposition of two known basis states, without making assumptions 
of the system-environment coupling. It is found that the obtained decoherence suppression can be 
highly sensitive to the ordering of the three UDD layers and can be remarkably effective with the 
correct ordering. The detailed theoretical results are useful for general understanding of the nature 
of controlled quantum dynamics under nested UDD. As an extension of our three-layer UDD, it is 
finally pointed out that a completely unknown two-qubit state can be protected by nesting four layers 
of UDD sequences. This work indicates that when UDD is applicable (e.g., when environment has 
a sharp frequency cut-off and when control pulses can be taken as instantaneous pulses), dynamical 
decoupling using nested UDD sequences is a powerful approach for entanglement protection. 

PACS numbers: 03.67.Pp, 03.65.Yz, 07.05.Dz, 33.25.+k 



I. INTRODUCTION 

Virtually all quantum systems are coupled to an envi- 
ronment and hence suffer from decoherence. Even more 
troublesome, the implication of decoherence for entan- 
gled states is far more severe than for a single quan- 
tum system. For example, even though a superposi- 
tion state of one qubit cannot be completely decohered 
within any finite time, the entanglement between two 
such qubits may be totally destroyed by decoherence 
within a very short time [1-3]. Evidently then, de- 
veloping useful schemes to protect quantum entangle- 
ment from environment-induced decoherence is crucial 
for entanglement-based quantum technologies. 

One promising approach towards decoherence suppres- 
sion is dynamical decoupling (DD) [4], which advocates 
the application of a sequence of instantaneous control 
pulses to effectively average out the system-environment 
coupling. Important extensions of the original DD ap- 
proach have also been developed, e.g., "concatenated dy- 
namical decoupling" pulses [5| , soft but optimized pulses 
under an energy cost constraint Q or a minimal leak- 
age requirement Q, and Uhrig's dynamical decoupling 
(UDD) [HH sequence that can achieve a very high effi- 
ciency of decoherence suppression, i.e., decoherence sup- 
pression to the iVth order with a sequence of N instan- 
taneous control pulses. Although UDD requires a sharp 
frequency cut-off in the bath spectrum plT[l3j and can- 
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not operate well if we set limitations on the control pulse 
width [l4| , UDD has attracted substantial interests soon 
after its discovery. In addition to its high-order suppres- 
sion of decoherence (at least in theory), UDD for single- 
qubit decoherence suppression could be powerful because 
it works for most general system-bath coupling 10], for a 
bath that has unknown spectral density (but with a sharp 
cut-off), and for time-dependent system-bath Hamiltoni- 
ans as well [15J . Experimental studies of UDD under two 
specific situations have been reported [l~6l - fl8j . For a re- 
cent concise review on UDD-related theoretical studies, 
seeRef. [H. 

So far the majority of DD studies have focused on 
single-qubit systems [2Cj . Hence it is urgent to investi- 
gate if quantum entanglement, e.g., two-qubit entangled 
states, can be well protected by DD. Because preserving 
two-qubit entanglement is more subtle than preserving 
single-qubit coherence, we wish to find a general control 
scheme to achieve good entanglement preservation with- 
out assuming any particular form of system-environment 
coupling. Indeed, given so many different ways of cou- 
pling a composite quantum system to an environment, in 
general a specific assumption about system-environment 
coupling may over-simplify the issue of entanglement pro- 
tection. As such, a universal and efficient DD scheme 
for entanglement protection should be of sufficient in- 
terest. Certainly, if under certain environments some 
crucial information about system-environment coupling 
becomes available, then a general DD scheme may be 
further reduced, a situation exploited in the first exper- 
imental study of entan glem ent protection by DD in a 
solid-state environment 1211] . 
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In our early study 22] , it was shown that an extended 
UDD can also lock a known two-qubit entangled state to 
the iVth order with N control pulses. The present study 
is concerned with the protection of unknown entangled 
states, a situation that is more relevant for quantum in- 
formation processing. To that end we extend the recent 
work by West et al. [23j], where a scheme based on two 
layers of UDD sequences is proposed to suppress both 
the population relaxation and transverse dephasing of 
one single qubit. In particular, we show that by nest- 
ing three layers of UDD sequences, it is possible to lock 
a two-qubit entangled state as an unknown superposi- 
tion of two basis states, to the Nth. order, using about 
A 3 control pulses in total. This entanglement protec- 
tion scheme is independent of how the two-qubit sys- 
tem is coupled to its environment. The control operators 
used in our nested UDD are also explicitly constructed 
based on two arbitrary, but known, basis states. We 
shall also show that there are different scenarios in con- 
structing the control operators. Interestingly, it is found 
that the ordering of the nested UDD layers can be a 
crucial factor for achieving high-order entanglement pro- 
tection. This intriguing ordering dependence is absent in 
two-layer UDD for single-qubit decoherence control, thus 
offering more insights into quantum decoherence control 
via nested UDD. 

As a further extension of our three-layer UDD, we show 
that it is possible to construct a four-layer UDD scheme 
such that a totally unknown two-qubit state can be pro- 
tected using about A 4 pulses of single-qubit control op- 
erators. Experimentally, this might be even more chal- 
lenging than realizing A 3 pulses of two-qubit operators. 
However, such a theoretical possibility is hoped to moti- 
vate future studies. 

This paper is organized as follows. In Sec. II, after in- 
troducing the most general system-environment coupling 
for two-qubit systems, we briefly review our previous ex- 
tension of UDD from one-qubit to two-qubit systems. 
Emphasis is placed on the key requirements to achieve 
such an extension. In Sec. Ill, we consider two schemes 
for nesting UDD in three layers in order to protect un- 
known two-qubit entangled states. Supporting numerical 
results are also presented. Section IV discusses a four- 
layer UDD scheme, followed by Sec. V that concludes 
this paper. 
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For convenience each term in the above total Hamiltonian 
is assumed to be time independent (this assumption can 
be lifted). Here Co represents the self-Hamiltonian of the 
bath, a x , a y , and a{ are the standard Pauli matrices for 



the first (j = 1) or the second (j = 2) qubit, 



U 7J 



and 



Cjs (7, <5 = x, y, z) represent arbitrary smooth bath oper- 
ators. From Eq. (1), it is seen that in general, the bath 
may interact with each individual qubit, or modulate the 
mutual interaction between the two qubits. The latter 
situation naturally arises if, for example, the bath can 
induce phonon excitations in a solid and hence perturb 
the relative distance between the two qubits embeded in 
the solid. The above total Hamiltonian is in the most 
general form, because it can be regarded as a linear ex- 
pansion over all possible 16 linearly independent basis 
operators operating on a four-dimensional Hilbert space, 
with the expansion coefficients containing arbitrary bath 
operators. Note however, the frequency spectrum of the 
bath is assumed to have a hard cutoff so that the general 
theory of UDD is applicable. 



II. PROTECTING A KNOWN TWO-QUBIT 
STATE BY UDD 

A. General total Hamiltonian of a two-qubit 
system interacting with a bath 

In terms of system-environment coupling, two-qubit 
systems are far more complex than one-qubit systems. 
A general total Hamiltonian describing a two-qubit sys- 



The basis operators used in the above-mentioned ex- 
pansion can be taken as {-Rj}i=i,2,--- ,16 = Wk ® c;}, 
with <7fc, <ji € {/, o x , &y , &z } {I the unity operator for the 
two-qubit Hilbert space) and the orthogonality condition 
Tr(RjRk) = 45 jk- This choice of basis operators is rather 
arbitrary. Purely for the sake of discussions below, we 
find it convenient to define two new sets of basis oper- 
ators. Let |0), |1), |2), and |3) be the four orthogonal 
basis states of the two-qubit Hilbert space, we define the 
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following two new sets of basis operators, 
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r ls =yi 5 = |o>(i| + |i><o|, 

Y 16 =Y 16 = -i(\l){0\-\0){l\). 

Then our general total Hamiltonian H may be re- 
expressed as 

16 16 

H = Y,WiYi=Y,W k Y k , (3) 

i=l fc=l 

where Wi or W k are the associated new expansion co- 
efficients containing bath operators. The motivation of 
using Yi (i = 1 — 16) as the basis operators is purely for 
convenience. For example, most of them have only one 
nonzero matrix element and some others are chosen to 
form a standard SU(2) subalgebra. Similarly, the basis 
operators Yi (i = 7 — 14) are chosen to simplify our later 
calculations involving states (|0) + |l))/\/2 [see Eq. (fj"9|) ]. 
Note also that the new basis operators defined in Eq. ([2]) , 
which still form the generating algebra of H, may not 
take a Hermitian form. This is not an issue because their 
linear superpositions still generate all possible Hermitian 
operators for a two-qubit system. 



B. Locking a known two-qubit state by an 
extended UDD scheme 

Given a known but arbitrary two-qubit state, here as- 
sumed to be |0) without loss of generality, we first con- 
struct a control operator 

X o = 2|0)(0|-J, (4) 

with Xq = I. As recently pointed out in Ref. 19], such 
a control operator was also considered in Ref. [2J| before 



UDD was discovered. We can now split H into two parts, 
H = i?o + Hi , 

10 

H = J2 W ^ 

i=l 
16 

Hi = ^WiYi, (5) 

i=ll 

with the commuting relation [Xq, Hq] = 0, and the anti- 
commuting relation {Xq,Hi} + = 0. We proposed in 
Ref. [lH the following control Hamiltonian describing a 
sequence of extended UDD 7r-pulses over a duration of 
T, i.e., 

N X 
H^^Sit-T^, (6) 

with the UDD timing Tj given by 

T 7 - =Tsin 2 (— ^— ), j = 1,2- •■ ,N. (7) 

For odd TV, an additional control pulse is applied in the 
end. Then the unitary evolution operator for the whole 
system of two qubits in a bath for the period t = to 
t = T is given by (K = 1 throughout) 

U N (T) = X^e-^ H ° +H ^ T - TN \-iX Q ) 

x e -^ffo+H 1 ](T N -T N _ 1 )(_ zXo ) 

x e -iWo+m]{T 3 -T 2){ _ iXo) 

Exploiting [ff ,-V ] = and {Hi,X } + = 0, one can di- 
rectly use the UDD universality proof developed by Yang 
and Liu [ill Oil , yielding 

U N {T) = U c N vcn + 0(T N+1 ), (9) 

where 

+oo 

U c N vcn = cxpH^oT) ^(-z) 2fc A 2fe , (10) 

k=0 

with A 2 k only containing even powers of H[(t), defined 
by H[ (t) = exp(iH t)Hi exp(-iH t). 

Because {Hi,X } + = 0, one has {H{(t), X } + = 0. 
As such, any even power of H\ (t) will commute with 
X , e.g, [H((ti)H[(t 2 ) 7 X } = H[(t 1 ){H((t 2 ).X }+ - 
{H[(ti), Xo} + H[(t 2 ) — 0. This important observation 
indicates that A 2k can be expanded as a linear superpo- 
sition of all possible basis operators that commute with 
X . That is, 

10 

A 2k = J2 A i y h (11) 

i=l 
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where Ai are the expansion coefficients containing bath 
operators. Clearly then, to the Nth. order, Ujq{T) can 
be expressed as a combination of Yi, Y2, • • • , Y w only. 
Using the closure of this set of operators, i.e., 

(10 \ / 10 \ 10 

^AiYA x>r fe = ( 12 ) 
»=i / \fc=i / (=1 

we further obtain 

U N {T) = exp(~ lj ff c u ff DD - 1 T) + 0(T N+1 ), (13) 

where 

10 
i=l 

with D\ i being the expansion coefficients. 

The outcome of applying a UDD sequence of Xq is now 
evident by comparing the original total Hamiltonian H 
in the absence of control with the effective Hamiltonian 
i/™ 0-1 realized by UDD. In essence, the UDD sequence 
based on Xq efficiently removes the operators Y\\, Y12, 
• • • , Y\q from the initial generating algebra of H, thus 
suppressing all possible coupling between the pre-chosen 
state |0) and all other states. Therefore, given a known 
entangled state |0), we can protect this state to the Nth 
order with N [or (N + 1)] instantaneous control pulses, 
a result analogous to single-qubit UDD. 

For later use, we list below three key requirements in 
achieving a UDD-reduced effective Hamiltonian H^f ' 1 
to the Nth order, from a general Hamiltonian describing 
two qubits plus a bath: 

(i) Construction of a control operator (e.g., Xq) whose 
square equals the unity operator. This control op- 
erator will be used to form a UDD sequence of N 
instantaneous pulses [e.g., Eq. ©]. 

(ii) Separation of the bare system-bath Hamiltonian 
into two terms, say Hq and H%, with Hq commuting 
with the control operator and Hi anti-commuting 
with the control operator. 

(iii) Algebra closure of the operators forming Hq , which 
becomes the generating algebra of a UDD-reduced 
effective Hamiltonian. 

In our following considerations we will make a number of 
references to these three requirements. 



III. PROTECTING UNKNOWN TWO-QUBIT 
ENTANGLED STATES BY NESTED UDD 

A. First scheme for nesting three UDD layers 

The Xq control operator in Sec. II is based on the 
knowledge that the state to be preserved is |0). It is 



even more useful if we can develop a scheme to protect 
unknown two-qubit entangled states. To that end, let 
us assume that an unknown two-qubit state |^(0)) is a 
superposition of two orthogonal basis states, e.g., |0) and 
|1). Though our considerations below are general, to be 
specific we assume |0) = | ft)) an( l |1) = I 44) 1 where 
t and 4- represent spin-up and spin-down states of each 
qubit. The unknown two-qubit state to be protected can 
be written as 

|V(0)) = a\0)+p\l) 

= a|tt) + /3|W), (15) 

where the two unknown coefficients a and (3 satisfy |a| 2 + 
|/3| 2 = 1 at time zero. Can we efficiently protect such type 
of unknown entangled states by further extending UDD? 
Note that this problem is different from a single-qubit 
case because the population can leak out from the initial 
two-dimensional subspace. 

We use p(t) to represent the total density matrix of 
the system and the bath at time t, evolving from a di- 
rect product state of IV'(O)) and some initial state of the 
bath. The protection of the state [^(O)) requires to freeze 
multiple coherence properties, i.e., diagonal populations 
Tv[p(t)\ ttXtt |] « M 2 and Tr[p(i)| U)(U |] a |/3| 2 , 
as well as the off-diagonal phase property Tr[p(i)| ft) (44 
|] sa a/3* . This motivates us to extend the UDD nest- 
ing scheme in Ref. [23[, where both single-qubit popu- 
lation relaxation and single-qubit transverse dephasing 
are suppressed in a near-optimal fashion. Certainly, our 
problem here is more demanding: in single-qubit sys- 
tems with a two-dimensional Hilbert space, the locking 
of one projection probability onto one basis state auto- 
matically freezes the projection probability onto a sec- 
ond basis state, whereas here the locking of the diago- 
nal probabilities at |a| 2 and |/3| 2 should be respectively 
achieved by control pulses. Given that two layers of UDD 
sequences are needed for complete single-qubit decoher- 
ence control in Ref. 23] , it is a natural guess that we will 
at least need three layers of nested UDD sequences. 

We can now directly make use of our results in the 
previous section. In the first step, we lock the diagonal 
property Tr[p(t) \ ttXtt I]- This can be achieved by con- 
sidering an innermost layer of UDD sequence of Xq , such 
that all possible coupling between |0) and all other or- 
thogonal states can be efficiently removed. Doing so, we 
reduce H to i?™ 0-1 to the iVth order, as elaborated in 
Sec. II. The population on state |0) is hence locked. 

In the second step, we treat a decoherence control 
problem for the effective Hamiltonian We as- 

sume that i/™ 0-1 resulting from the innermost layer of 
UDD is a sufficiently smooth function of time, such that 
a second layer of UDD can be applied to H^f®' 1 . Note 
however, though this is a natural assumption and Ref. 

reasoned about the smoothness of analogous UDD- 
reduced effective Hamiltonians, some counter examples 
might exist fl9j . With this smoothness assumption exer- 
cised with caution, we now introduce a second UDD layer 
to lock the second diagonal property Tr[p(t)| 44) (44 |]. 
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As is clear from Sec. II, in order to remove all possible 
couplings between |1) and all other states, one may ap- 
ply the control operator X% = 2|1)(1| — / with X 2 = /. 
To examine if this is feasible, we decompose H^f®' 1 into 
two terms, i.e., 

10 

2 = 1 

_ rrUDD-1 , ttUDD-1 

with 

6 

»=i 

10 

i=7 

Because [H^- 1 , X{\ = and {H^~ 1 ,X 1 }+ = 0, it 
is seen that H^f ' 1 and X\ guarantee requirement (ii) 
outlined in Sec. II. Finally, it is straightforward to see 
that the operators Yi, i = 1 — 6 form a closed algebra [re- 
quirement (iii) above] . We hence expect that if a second 
layer of UDD sequence of Xi is applied, to the Nth. order 
the dynamics of H^ f P B ~ 1 under the second control layer 
becomes that of a simpler effective Hamiltonian (denoted 
H^f B ~ 2 ) generated by a further reduced algebra. That 
is, 

H™»- 2 =Y.D 2 ^ (18) 

i=l 

where D^^i represent the expansion coefficients due to 
two UDD layers. 

Examining the self-closed set of operators that form 
H^f ' 2 , one sees that only the component -02,6^6 = 
£*2,6[|0)(0| — can affect the intial superposition 

state a\0) + j3\l). Further, this Yq component does not 
change the populations on states |0) and |1), so it rep- 
resents a pure dephasing mechanism. To efficiently sup- 
press this pure dephasing, we now consider a third, out- 
ermost UDD layer. Assuming again that iJ^ 00 " 2 is suf- 
ficiently smooth for UDD to apply, we introduce the fol- 
lowing "phase" control operator 

X, = [|Q) + |1)][<0| + <1|]-J, (19) 

with X^ = I. Separating D " 2 into two terms, we 
obtain 

rrUDD-2 rrUDD-2 , ttUDD-2 / on \ 
H eS - H e S,0 + H cS,l > 

with 

5 

#rff D o D ~ 2 = J] D2,jYj] 
i=i 

H™»- 2 = D 2)6 Y 6 . (21) 



Interestingly, [iJ c u ff D D " 2 , X ] = and {ff™ " 2 , A }+ = 
0. Further, the operators Yi, i = 1 — 5 form a closed 
algebra. All the three key requirements for UDD are 
again met for this outermost layer. The final reduced 
effective Hamiltonian after three layers of UDD is hence 
formed by five operators, i.e., 

i? c u ff DD - 3 -E^, (22) 
»=i 

where D^ i represent the expansion coefficients due to 
three UDD layers. Referring to the subspace spanned by 
|0) and |1), H"JgP D ~ 3 only contains an identity operator for 
that subspace. Hence any unknown initial superposition 
state a|0) + /3|1) is well preserved to the iVth order. 

In terms of step- by-step reduction of the generating al- 
gebra associated with the effective Hamiltonians at each 
level, the following chart summarizes how our nesting 
scheme reduces a general total Hamiltonian to a much 
simplified form: 



Yi, i = 1,2,- ■■ ,16 
JjJfo, UDD-1 

Yj, i = l,2,--~W 
JJJfi,UDD-2 

Y u z = l,2,--~6~ 
4A ,UDD-3 

Yi, i = l,2,--- ,5 . 



The operators inside each box represent the elements of 
the generating algebra before or after a certain layer of 
UDD. 

For completeness, we also explicitly present here the 
timing of the control pulses within each UDD layer. In 
particular, the control operator Xq in the innermost layer 
is applied at 

T j<k ,i = T jtk + (T.j^k+i - T jtk ) sin 2 [^p^j , (23) 

where Tj^ is the UDD timing for X\ in the middle layer. 
Tj >k is given by 

T jtk = Tj + (T j+1 Tj) sin 2 (^g) > ^ 

where Tj represents the timing for A^ in the outermost 
layer and it is already given by Eq. ([7]). Similar to single- 
qubit cases, for each layer, if N is odd then an additional 
control operator is applied at the end of each sequence. 
Overall, N 3 [or (N + l) 3 ] control pulses are applied to 
achieve decoherence suppression to the Nth order. Be- 
cause states |0) and |1) play a similar role here, our anal- 
ysis above equally applies if X\ and Xq are exchanged. 
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TABLE I: Explicit construction of control operators for 
nesting three layers of UDD sequences in order to protect 
unknown two-qubit entangled states q|0) + with 
|0> = | ft) and |1> = | U). 



where are the expansion coefficients. To consider 
the second UDD layer, we rewrite i/^ 00 " 1 in terms of 



Constructed Operator 


Explicit Form 
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{o\o'l + al +ol-l) /2 


Xr 


(aWl - al - 0* - I) /2 


X^ 


{(t\o% — (j\a'l + a\a 2 z - I) /2 


Xqa 





To realize such a UDD nesting scheme with this high- 
order decoherence suppression, the involved two-qubit 
control operators are nonlocal control operators in gen- 
eral. Taking |0) = | tt) an d |1) = | 11) as an example, 
Table I lists the three control operators (another operator 
to be explained later) in terms of the familiar Pauli ma- 
trices. The identity operator in the expressions for Xq, 
X\ and Xj, in Table I is not important. Experimentally, 
realizing such two-qubit control operators is analogous 
to realizing quantum computation in a two-qubit system. 
The true challenge might lie in realizing a sufficient speed 
of such two-qubit operations. 



B. Wrong ordering of three UDD layers 

At this point an interesting question arises. That is, 
does the ordering of the three nested UDD sequences 
matter or not? To answer this question let us first ex- 
change the ordering of the two sequences of Xq and 
X<p, such that X^ is in the innermost layer and Xq 
is in the outermost layer. We denote this ordering as 
Xq — X\ — X$. In the following we shall stick to this 
convention for ordering, i.e., an operator appearing at 
the rightmost (leftmost) will be placed in the innermost 
(outermost) layer. Because the X^ layer is now operat- 
ing directly on the bare Hamiltonian H, we re-partition 
H in Eq. ([3]) into the following two terms, i.e., 

/ 5 10 \ 



\k=l k=7 I 

+ (w e Y 6 + W 16 Y 16 + W k Y k \ , (25) 

V k=n ) 

where operators Y k are defined in Eq. ([2]) . The first term 
in the above equation commutes with X^, whereas the 
second term anti-commutes with X^. All the operators 
contained in the first term form a closed algebra. The 
innermost UDD sequence of X$ hences yields an effective 
Hamiltonian 

5 10 



fVUDD-1 



{y^D'^Yi + D'^Y^ , (27) 



where D[ i is connected with D\ k via a simple relation 
between Yi and Yi [see Eq. (J2J]. 

The middle layer of UDD sequence would be based 
on the control operator X\, Among those basis opera- 
tors that form _ff™ D_1 , X\ commutes with those in the 
first line of Eq. (|2"T|) and anti-commutes with those in 
the second line of Eq. ([2T]) . This indicates that the final 
evolution operator associated with this UDD layer can 
be cast into a form similar to Eq. (|10l) . We next investi- 
gate if the key requirement (iii) of UDD can be satisfied. 
Interestingly and somewhat unexpectedly, this is not the 
case: operators Yi with i = 1, 2, • • • ,5 together with Yi 
with i = 11 — 14 cannot form a closed algebra. For exam- 
ple, 2YnYi2 — Y 2 will yield Y 6 , which is already outside 
this collection of basis operators. As a result, the appli- 
cation of this second control layer will not yield a further 
reduced effective Hamiltonian. Such a nesting scheme 
then breaks down due to its wrong ordering! Indeed, in 
the correctly ordered case, the dephasing operator Y 6 will 
be suppressed by but here it resurfaces (from even 
powers of those operators that commute with X\) after 
the Xrf, layer is already applied. 

Note however, due to this wrong ordering, the unde- 
sired operators (such as Yq that cannot be suppressed by 
the outermost Xo layer) reemerge from multiplications of 
a set of basis operators. We hence intuitively expect its 
impact on decoherence control to be at least a second- 
order effect. We will come back to this when discussing 
our numerical results. 

One may wonder what happens to other ordering? Be- 
cause states |0) and |1) play the same role here, there is 
only one non-equivalent ordering left, i.e., X\ — X^ — Xq 
(or equivalently, X — X^ — X{). The effective Hamilto- 
nian after the innermost layer of UDD sequence is hence 
still given by H^fV-i = Y™ = i D ^ Y i- Since the next 
layer of control is X^, one checks if H^ l f f' u ~ 1 can be de- 
composed into two terms in accord with the properties 
of X,p [to fulfill requirement (ii) outlined above]. Inter- 
estingly, though this procedure can be easily done in the 
entire two-qubit operator space, it cannot be done here 
within the reduced generating algebra of i?^ 00 " 1 . That 
is, among the set of basis operators that form H^g >u ~ 1 , 
some operators or their arbitrary combinations are nei- 
ther commuting nor anti-commuting with X^. For ex- 
ample, 



i=7 



(26) 



[r 7 ,A^0; {Y 7 ,X^}+^0. 



(28) 
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So even without checking if there is a self-closed set of 
operators that commute with X^, it is already seen that 
the nesting scheme breaks down directly in this order- 
ing. Thus, out of three non-equivalent ordering of the 
three control operators considered here, only the order- 
ing advocated in the previous subsection may achieve 
high-order protection of two-qubit states. 

From a more general perspective, the dependence of 
the controlled dynamics on the ordering of the three UDD 
layers can be explained as follows. Consider a particular 
time interval AT for a time-independent Hamiltonian H , 
during which each of the two control operators A and B 
(A 2 = 1 and B 2 = 1) are applied twice with a certain 
ordering. If B is nested inside, the associated unitary 
evolution is given by 

U A - B = ABe-* HAT BA = e -i[(AB)H(BA)]AT. (2g) 

whereas for the other ordering, the unitary evolution is 
given by 

U B - A = BAe- lHAT AB = e -<K*A)ir(AB)]AT (30) 

As seen from the right hand side of Eqs. (f2U)l and 
the two ordering leads to two effective Hamiltonians 
(AB)H(BA) and (BA)H(AB). In general these two ef- 
fective Hamiltonians arc different, thus giving rise to the 
ordering dependence. Note however, if for a concerned 
Hilbert subspace AB = BA or AB = -BA, i.e., if the 
two control operators commute or anti-commute, then 
these two effective Hamiltonians are identical and hence 
the ordering dependence no longer exists. This further 
explains why we cannot exchange the ordering between 
and Xq or the ordering between X^ and X\ , but can 
exchange the ordering between Xq and X\ . 

C. Alternative Nesting Scheme 

For the two-dimensional subspace spanned by states 
|0) and |1), any two orthogonal states |0') and |1') can 
be adopted to construct analogous control operators Xqi , 
Xy, and X^i for three UDD layers. Therefore, when 
it comes to an actual implementation, there are infinite 
possibilities to realize three nested UDD layers in order 
to protect an unknown entangled state in that subspace. 

Even for fixed basis states |0) and |1), the nesting 
scheme with the right ordering analyzed above is not the 
only solution. Consider the following control operator 

X ,i = 2[|0)(0| + |1)(1|] - I = X + Xx + I (31) 

with Xq i=I- The explicit form of Xo,i is also given in 
Table I. Following the analysis in Sec. II, it is straight- 
forward to confirm that a UDD sequence of A^i can effi- 
ciently freeze the total population in the two-dimensional 
subspace spanned by |0) and |1). We can now construct 
an alternative nesting scheme using Xq i , X^ and X\ (or 

*o). 



The ordering of X^ — X\ — Xq,i is studied first. It 
is found that each of the three UDD layers satisfies the 
three requirement outlined in Sec. II and hence yields a 
simple effective Hamiltonian to the Nth order. Qualita- 
tively, it is also obvious why this scheme is expected to 
work. The innermost layer effectively lock the population 
in a two-dimenional Hilbert subspace. Then, in essence, 
the next two layers are similar to those in single-qubit 
two- layer UDD |23], insofar as a second layer locks one 
population and the outermost layer freezes the relative 
phase between the two projection amplitudes. Extend- 
ing this analogy, it is expected that the ordering of the 
middle layer and the outermost layer can be exchanged 
(indeed, the associated control operators anti-commute 
in the two-dimensional subspace). This can be more for- 
mally analyzed by working out the detailed algebra layer- 
by- layer. In particular, for the ordering of X^—Xi—Xo^, 
the generating algebra for the effective Hamiltonians af- 
ter each control layer reduces in the following fashion: 



Yj,i = l,2,--- ,16 
tyX ,i, UDD-1 
Yj, i = l,2,---,6;Yi 5 , Yj 6 
JJJfi,UDD-2 
Y u i = l,2,--- ,6 

^A ,UDD-3 
Y%, i = 1, 2, • • • ,5 . 



For the ordering of X\ — X^ — Xq,i, we have 



Yj,i = l,2,--- ,16 
JJJf ,i, UDD-1 
Y,i = l,2,--- ,6; Yi 5 , Yj 6 
4X ,UDD-2 
Yj,i = l,2,--- ,5;F 15 
JJJfi,UDD-3 
Yi, i = 1,2,- ■■ ,5 . 



As seen above, both nesting strategies are successful and 
in the end the same generating algebra for the final ef- 
fective Hamiltonian is reached. 

Can we place a sequence of Xo t i in the middle layer 
instead? For the ordering of X^ — _X"o,i — X\, since the 
population on state |1) is already locked by the inner 
most layer, the role of the second layer is equivalent to 
further locking the population on state |0), thus playing a 
similar role as Xq. The outermost layer then freezes the 
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relative phase between the two projection amplitude 
analogous to the case of — X — X x . In this sense, tl 
ordering of X^ — X0.1 — X\ does not provide anythir 
new. Indeed, based on the discussion at the end of tl 
last subsection, because [JCo,i,Xi] = 0, their ordering 
expected to be exchangeable. 

Consider then the other ordering X\— Xo,i — X^ (whk 
is also equivalent to Xq — X^i — X^). In this case, tl 
effective Hamiltonian reduced by the innermost layer 
given by H^P ' 1 in E T 03- To analyze the effect of tl 
second layer, we re-split -ff™ 13 " 1 in Eq. (l27l) as 



fVUDD-l 
-"eff 



\J2D' lti Y i + D' 1A5 Y 15 j 

14 



io- 6 



i=7 



I 1 I ■ t 1 I 1 



X -X -X„ 

it 1 o 




where the first term commutes with Xo,i and the secon„ 
term anti-commutes with A"o.i- As an interesting out- 
come, now the set of operators that commute with the 
second-layer control operator also form a closed algebra, 
thus paving the way for the third UDD layer. Indeed, 
the outermost Xo layer further reduces the algebra by 
removing the Y15 component and hence yields an effec- 
tive Hamiltonian seen before. Summarizing, the explicit 
algebra reduction route for X\ — Xq^i — X^ becomes 



Yi, i = l,V.. ,16 



11 A^, UDD-1 



Yi, i = !,-■■ ,5; Yi, i = 7,8,- ■■ ,15 



^X ,i,UDD-2 



Y u i= 1,2, •-• ,5;Y, 



l.r, 



J^Xi,UDD-3 



Yi, i — 1, 2, • • • , 5 



We have also examined what happens if the Xq.i layer 
is placed in the outermost layer. For reasons analogous 
to our first nesting scheme using Xq, X\ and X$, such 
type of ordering cannot simplify the effective Hamiltoni- 
ans layer-by-layer. This concludes this subsection. 



D. Numerical study 

Similar to our previous work (22j |. we use a five-spin 
system to carry out simple numerical experiments. Two 
of the five spins are identified as our two-qubit system, 
and the other three spins are regarded as the bath. To 
avoid assumptions about how the system is coupled to the 
bath, we work with the following general total Hamilto- 



FIG. 1: (color online) Averaged trace distance D between ini- 
tial and final two-qubit states in the presence of three-layer 
nested UDD. TV is the number of instantaneous control pulses 
for each layer. The three control operators Xo, Xi, and X$ 
are explicitly constructed in Table I. The total system-bath 
Hamiltonian is modeled by Eq. (|33f) . Note a strong depen- 
dence on the ordering of the three UDD layers. Top two 
curves are for Xo — X<j, — X\ and Xi — X^ — Xo ; bottom two 
curves are for X<p — Xq — X\ and — X\ — Xq. All the vari- 
ables plotted here and in all other figures are in dimensionless 
units. 



H = E E 

5 5 

+ E E E E c T<°s ( 33 ) 

m=l j={x,y,z} n>m 6={x,y,z} 



in dimensionless units, where all the coefficients 6 7>m 



and 



c^g 1 are randomly sampled from the range [—0.5,0.5]. 
We average our results over ten random realizations of 
this five-spin system-bath Hamiltonian. In addition, to 
demonstrate that our approach does not depend on the 
actual form of an initial superposition state a\ ft) +0\ 44- 
), we further average our results over ten initial states 
with randomly sampled coefficients a and /3 under the 
constraint \a\ 2 + \/3\ 2 = 1. 

Figure 1 depicts the averaged trace distance (denoted 
D) between the system's reduced density matrix at time 
t = 0.1 and its initial state, for N = 1 — 10. Different 
ordering of three UDD sequences X , X\, and X^ are 
plotted together for comparison. Consider first the two 
cases (bottom two curves) with the correct ordering, i.e., 
X^ — X\ — Xq and X^ — Xq — X\ . For these two cases a 
remarkably high fidelity is achieved in locking the initial 
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unknown superposition state. For N = 10 (totally N 3 
UDD pulses), D already reaches the 10 -10 level. The al- 
most linear scaling of log(£>) vs N is consistent with the 
expectation that the extent of the decoherence suppres- 
sion for a working nesting scheme is to the iVth order. 

Turning to the top two flat curves associated with the 
wrong ordering Xi — — Xq and X — X^ — X\ . Their 
D values do not decrease with N and stay at about 10 _1 . 
Therefore, for these two cases the three-layer nested UDD 
does not work at all due to the wrong ordering. This di- 
rectly confirms our early insights into the issue. In par- 
ticular, from the algebra considerations we observe that 
for the incorrect ordering here, the second layer of UDD 
directly breaks down because the effective Hamiltonian 
reduced from the inner most layer does not meet require- 
ment (ii). 

Still referring to Fig. 1, let us now discuss the middle 
curves associated with another type of wrong ordering, 
i.e., X\ — Xq — Xcf, and X — X\ — X$. It is seen that 
their D values first decrease and then tend to saturate 
as N increases. The smallest D values for N = 10 is 
about 10 -4 , which is about six orders of magnitude larger 
than in previous correctly ordered cases. However, this 
performance is at the same time better than the top two 
flat curves. As such, the wrong ordering here represents a 
weak deviation from an ideal nesting. This is consistent 
with our early intuition that for the current ordering, 
the UDD nesting scheme breaks down due to a high- 
order effect, i.e., the non-closure of a set of operators 
that commute with the control operator in the middle 
layer. 

As a comparison with a conventional dynamical decou- 
pling approach based on control pulses equally spaced 
in time, Fig. 2 presents the parallel results if, within 
each layer, the control operator is applied periodically. 
Clearly, in this case, irrespective of the ordering of the 
control operators, the performance of decoherence con- 
trol for N = 10 in all cases is about nine orders of mag- 
nitude worse than the best two cases in Fig. 1. It is 
also observed that the D values are very weakly depen- 
dent on N. Results here remind us that in addition to 
the ordering of the three layers, the timing of the control 
operators is essential. 

Finally, results for an alternative UDD nesting scheme 
based on AT ,i, X\, and X^ are shown in Fig. 3. The top 
two curves are for two incorrect ordering AT ,i — X\ — X^ 
and Xq^i — X^ — X±, with their D values saturating 
at about 10 -4 as N increases. It is noted that the 
A^i — X\ — X^ curve is similar to the Xq — X\ — X^ 
case in Fig. 1. This is understandable because the un- 
derlying mechanism for unsuccessful nesting is the same. 
However, the Xo,i — X$ — X\ case here has better per- 
formance than the Xq — X$ — X\ case in Fig. 1. This 
is interesting because in both cases, the effective Hamil- 
tonian reduced from the innermost layer does not satisfy 
requirement (ii) for the second layer. Clearly then, for 
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FIG. 2: (color online) Same as in Fig. 1, except that the con- 
trol pulses are applied periodically within each control layer. 
Note that the scale of the plotted D values is many orders of 
magnitude different from that in Fig. 1. 

the outermost layer, a sequence of Xq^i turns out to be 
superior to a sequence of X . This is somewhat expected 
because A"o,i in the outermost layer can still freeze the to- 
tal population in the two-dimensional subspace whereas 
Xq cannot. 

All other four curves in Fig. 3 display a roughly linear 
scaling of log(£>) vs N, indicating the success of three- 
layer nested UDD. Indeed, the layer ordering associated 
with these curves is all predicted to be correct in our 
theoretical analysis above. Interestingly, for fixed N, 
the performance for the ordering of X^ — X\ — X i or 
X^ — Xq x — X\ can differ from that for Xi — X^ — X i or 
X\ — Ao.i — X$ by about two orders of magnitude. Com- 
paring with the best performance here with that in Fig. 
1, a difference about two orders of magnitude is also ob- 
served. These numerical details indicate that even with 
the same scaling with N, the actual performance of a 
correctly ordered three-layer UDD may depend on the 
specific algebra reduction route. More insights into this 
intriguing finding might help to further understand the 
nature of decoherence dynamics under nested multi-layer 
UDD. 



IV. FROM THREE-LAYER UDD TO 
FOUR-LAYER UDD 

We are optimistic that for N - 10 a total of iV 3 UDD 
control pulses as proposed in this work might be achiev- 
able in some systems in the near future. If this is the 
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FIG. 3: (color online) Same as in Fig. 1, except that the three 
control operators used in the three UDD layers are now Xo,i, 
Xi, and X$. Top two curves are for Xo,i — X\ — X$ and 
Xo,i — Xtj, — Xi; bottom two curves are for X\ — JQ, — Xo,i 
and Xi — Xqa — X<j,. 



case, then as the next step one wonders if there exist 
even better schemes, at least in theory. In particular, as 
a result of three correctly ordered UDD layers, the gener- 
ating algebra for the final effective Hamiltonian H^f®' 3 
contains only five operators. Can we construct better 
control operators to reduce the algebra more rapidly? 
Can we even consider one more UDD layer to remove all 
possible system-environment coupling? 

As a brief summary of our latest progress along these 
two questions, we first note that, the two nesting schemes 
proposed in Sec. Ill treat the subspace spanned by |0) 
and |1) differently than the subspace spanned by |2) and 
|3). Indeed, we have assumed that the initial state is a 
superposition of states |0) and |1). If the initial state 
is totally unknown, then it is best to construct control 
operators that are symmetric with respect to the two 
subspaces. Upon completion of our studies of the two 
nesting schemes discussed in Sec. Ill, we find that the 
following three symmetry-adapted control operators can 
form another nesting scheme for three- layer UDD, i.e., 

Z X = |0><0| + |1)<1|-|2)<2|-|3><3|=Xo, 1 

Z 2 = |0)(0|-|l)(l| + |2)(2|-|3)(3| 

Z 3 = |0)(l| + |l)(0| + |2)(3| + |3)(2|. (34) 

That is, for each UDD layer, the three key requirements 
of UDD outlined in Sec. II are satisfied. Dramatically, 
the form of the three control operators Z%, Z 2 , and Z% 
is explicitly symmetric with respect to an exchange be- 
tween the (|0), |1)) subspace and the (|2), |3)) subspace. 



So what happens to the first subspace also applies to the 
second subspace. Their physical meaning is also clear: Z\ 
locks the total population within each subspace spanned 
by |0) and |1) or by |2) and |3), Z 2 locks the individ- 
ual populations on each state, and Z 3 finally suppresses 
the pure dephasing within each of the two subspaces 
of dimensional two. Therefore, this symmetry-adapted 
three-layer UDD scheme should have more efficiency in 
reducing the algebra layer-by-layer. For the ordering of 
Z-& — Z 2 — Z\, the associated algebra reduction route is 
found to be: 



Yi, i = 1,2,- ■■ ,16 



-U-Zi, UDD-l 



Yi, * = 1,2,... ,6; Yis, Y 



16 



JJ.Z 2 ,UDD-2 



Y u Y 2 , Y 3 , Y 6 



J|Z 3 ,UDD-3 



*i, Y 2 . 



The final effective Hamiltonian H^®' 3 after such three 
UDD layers is hence a combination of only two oper- 
ators: Y 2 and the unity operator Y\ (or equivalently, 
|2)(2| + |3)(3| and Y x ). Note that the dephasing be- 
tween the (|0), |1)) subspace and the (|2), |3)) subspace 
is still not suppressed. Perhaps even more remarkable, 
for these three symmetry-adapted control operators, they 
either commute or anti-commute, and consequently dif- 
ferent orderings of Z±, Z 2 and Z% can produce the same 
final generating algebra. 

One can further rewrite H^®' 3 in a more enlightening 
and symmetric form, i.e., 



rrUDD-3 
n cS;Z 



D 



3.2 



|1)(1|-|2)(2|-|3)(3|].(35) 



This finally brings us to our last theoretical question: 
can we further reduce H^^~ 3 by adding one more UDD 
layer? Our answer is yes in theory. This is somewhat 
obvious if one introduces the fourth control operator 



Z, ee |0)(2| + |2)(0| + |1)(3| + |3)(1| 



(36) 



with Z\ = 1. Clearly, because Z4 anti-commutes with 
the second operator in Eq. ([35]) and commutes with the 
identity operator, the fourth UDD layer based on Z4 
will finally yield an effective Hamiltonian as a certain 
bath operator multiplied by a unity-operator in the four- 
dimensional Hilbert space for a two-qubit system! So 
it is theoretically possible to protect a totally unknown 
two-qubit state against most general system-environment 
coupling, using about iV 4 control pulses in total. Further, 
since this four-layer scheme is intended to lock any two- 
qubit state, one may now arbitrarily choose the four or- 
thogonal basis states in order to simplify the four control 
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operators (this is not allowed in three-layer UDD because 
a two-dimensional subspace is chosen beforehand). In 
particular, if we consider a new set of basis states differ- 
ent than above, e.g., |0) = | tt),Jl) = | ft), |2) = Ut>, 
and 1 3) = | ft), then one obtains Z\ — a\, Z2 = erf, 
Z3 = erf, and Z4 = a\, which are only local control 
operators in this new representation. Such a four-layer 
solution is also numerically checked. 

V. CONCLUSION 

With both theoretical analysis and numerical study, 
we have shown how nested three-layer UDD can protect 
unknown two-qubit entangled states as a superposition 
of two known basis states, to the iVth order with about 
N 3 control pulses, without assuming a specific form of 
system-environment coupling. This is of much interest 
to current theoretical investigations of entanglement pro- 
tection. As a remarkable side result, it is found that the 
ordering of the three UDD layers can be a crucial factor. 
Numerical results support our theoretical considerations. 

Given the theoretical feasibility of extending UDD be- 
yond single-qubit systems, decoherence control via nested 
UDD should be of experimental interest as well. Though 
a rigorous mathematical foundation for nested UDD is 
still under development [III [25j], the success of three- 
layer UDD demonstrated here in two-qubit systems fur- 



ther strengthens the view that nested UDD can be a good 
strategy for decoherence suppression. Our approach can 
also be extended to protect an unknown superposition of 
two known basis states in an arbitrary multi-level system. 

As a final extension, in Sec. IV we also discussed how 
a totally unknown two-qubit entangle state can be pro- 
tected by applying four layers of UDD sequences that 
involve local operations only. Applying N pulses can be 
highly demanding in experiments, but the existence of 
such a theoretical solution should offer a useful reference 
point for future studies of entanglement protection. On 
the other hand, it becomes interesting to compare our 
main topic of this work, i.e., three- layer UDD, with the 
ultimate four-layer solution. The four-layer solution can 
be realized by local control operators but the required 
number of pulses may present experimental difficulties. 
By treating a particular class of two-qubit states, three- 
layer UDD can achieve entanglement protection of un- 
known states using much less control pulses, with the 
price that rapid nonlocal operations are required. 
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